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The characterization of the long-range order and fractal properties of DNA sequences has proved 
a difficult though rewarding task due mainly to the mosaic character of DNA consisting of many in- 
terwoven patches of various lengths with different nucleotide constitutions. We apply here a recently 
proposed generalization of the detrended fluctuation analysis method to show that the DNA walk 
construction, in which the DNA sequence is viewed as a time series, exhibits a monofractal structure 
regardless of the existence of local trends in the series. In addition, we point out that the monofractal 
structure of the DNA walks carries over to an apparently alternative graphical construction given 
by the projection of the DNA walk into the d spatial coordinates, termed DNA trails. In particular, 
we calculate the fractal dimension Dt of the DNA trails using a well-known result of fractal theory 
linking Dt to the Hurst exponent H of the corresponding DNA walk. Comparison with estimates 
obtained by the standard box-counting method allows the evaluation of both finite-length and local 
trends effects. 



I. INTRODUCTION 

The search for patterns in nature and their interpre- 
tation in terms of general principles is one of the main 
purposes of science. The explosive accumulation of DNA 
sequence data in the last two decades has provided a 
rich source of raw material that hides valuable hints 
about the evolutionary mechanisms of the genome or- 
ganization. Unveiling the patterns in those sequences 
has become an exciting challenge to the present gener- 
ation of statistical physicists and information scientists. 
In that vein, a truly remarkable result was the finding 
that intron-containing DNA coding regions exhibit long- 
range power-law correlations extending across more than 
10^ nucleotides, whereas intronless coding regions display 
only short-range correlations |^. 

Characterizing the long-range correlations in DNA se- 
quences is a highly nontrivial task because of the mo- 
saic structure of DNA consisting of patches with differ- 
ent nucleotide composition Q. In fact, in order to use 
the standard techniques (e.g., power spectrum and R/S 
fluctuation analysis) to study the correlations of DNA se- 
quences, it is necessary first to eliminate the local biases 
in nucleotide composition (trends), thus avoiding spuri- 
ous effects due to the mosaic character of the sequence. 
In the context of time series or records there are two 
main techniques to eliminate these trends, namely, the 
detrended fluctuation analysis (DFA) and the wavelet 
transform (WT) |4|. The former is a physically appealing 
ad hoc technique of easy implementation and the latter 
is a mathematically well-established transform used to 
study the regularity of arbitrary functions via the sys- 
tematic elimination of local polynomial behavior. Both 
techniques can be readily applied to the analysis of or- 
dered linear sequences, such as DNA, by considering the 
reading direction as the time axis (see e.g. |, | for the 
DFA and §, § for the WT analyses of DNA sequences). 

A common strategy to represent graphically a given 



DNA sequence consists of transforming it into a d- 
dimensional random walk, so-called DNA walk, by asso- 
ciating a space direction to each nucleotide (A,C,T,G) or 
class of nucleotides (e.g. purine and pyrimidine) and the 
time direction to the reading direction of the sequence 
Here, A, C, T and G are the bases adenine, cy- 
tosine, thymine, and guanine, respectively. In particu- 
lar, a very popular representation is the one-dimensional 
walk in which a purine (A or G) at position i is asso- 
ciated to one step down (jji — —1) while a pyrimidine 
(T or C) is associated to one step up {yi = 1). We 
note that the purine-pyrimidine classification is related 
to the hydrophobic-hydrophilic characteristics of the en- 
coded amino acids. Of course, there are many alternative 
DNA walks and the most complete one that considers 
the four bases equivalently in base space is the three- 
dimensional DNA walk in the sense that its projec- 
tions on appropriate axes or planes recover all possible 
one and two-dimensional walks. 

For the sake of concreteness, in this section we will 
focus on the fiuctuations of the cumulative variable of 
the purine-pyrimidine random walk only, defined as 



s(n) 



1, 



(1) 



where N is the length of the sequence. This variable 
is displayed in Fig. ^ for the intron-rich human beta- 
globin region (GenBank name HUMHBB) for which N = 
73326 as well as for a shuffled sequence with same length 
and nucleotide composition. The statistical quantity that 
characterizes the fluctuations of s{n) is the root mean 
square fluctuation F2{1) defined as 



F2(0= [As(/)]^- As(0 



1/2 



(2) 



where As{l) — s(n -\- 1) — .s(n) and the bars indicate an 
average over all positions n in the record. We expect F{1) 
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FIG. 1: Purine-pyrimidine random walk plot of the human 
beta-globin region (black line) and a shuffled sequence of the 
same nucleotide composition (grey line). 



to increase with increasing I, the length of the window 
considered. ExpHcitly, in the case of stationary series the 
fluctuations are described by a power law 



F2 (0 ~ I 



H 



(3) 



where the scaling exponent (Hurst exponent) H equals 
1/2 in the case of random and short-ranged correlated 
sequences. Any value different from 1/2 is evidence of the 
existence of long-range correlations in the sequence |l^, 
|lT| . In addition, for self-affine records (e.g. DNA walks) 
there is a simple relation between H and the local fractal 
dimension of the record, Dr = 2 — H. We note that the 
power law behavior implies that there is no characteristic 
length scale, i.e., very large fluctuations are likely due to 
the same mechanisms as smaller ones. 

Complex records are unlikely to be fully characterized 
by a single scaling exponent as evidenced by the case 
where the scaling behavior is different in distinct parts 
of the series so that H is actually dependent on the part 
of the record being analyzed. A complete framework to 
describe the situation where a multitude of scaling ex- 
ponents is required is provided by the multifractal for- 
malism ||ll|, |l^ and, in this sense, the multifractal spec- 
trum may be viewed as the ultimate tool to characterize 
a stationary time series. However the direct calculation 
of ^2(0 using Eq. (||) or the application of the standard 
multifractal formalism yield wrong results for nonstation- 
ary time series that are affected by local trends. In this 
contribution we apply a generalization of the DFA re- 
cently proposed by Kantelhardt et al. for the character- 
ization of nonstationary time series to study the multi- 
fractal properties of DNA walks jl^ . In agreement with 
a previous multifractality analysis based on the wavelet 
transform modulus maxima (WTMM) method Q, we 
find that the DNA walks exhibit a simple monofractal 
scaling behavior. The multifractality analysis of one- 
dimensional DNA-walks is the subject of Sec. ||. 




x(n) 



FIG. 2: Trail of the human beta-globin region (black line) and 
a shuffled sequence of the same nucleotide composition (grey 
line) using the two-dimensional representation T {xi = 1), A 
{xr = -1), G {y, = 1) and C [y^ = -1). 



A seemingly alternative graphical representation of 
DNA sequences that has received less attention and so 
less criticism than the aforementioned DNA walks is the 
so-called DNA pseudo random- walks 0, |l5[ ^ ^ 
They are simply the projection of the DNA walks into 
their d spatial components, i.e., the trails of the DNA 
walks. Of course, drawing these trails makes sense for 
d >2 only, and Fig. |^ illustrates such a trail for a d = 2 
representation of the HUMHBB sequence where right, 
left, up and down steps correspond to the presence of 
nucleotides T, A, G and C, respectively. This represen- 
tation preserves some basic symmetries of the DNA such 
as complementarity, reflection, substitution and compat- 
ibility [|l5[ . 

Despite the clear visual evidence of local trends in the 
trails exhibited in Fig. |^, it has been a common prac- 
tice to apply standard fractal and multifractal methods 
(e.g., box counting 1^ and sandbox |jl^) to character- 
ize their shapes without much concern about the mosaic 
structure of the underlying DNA sequence. Actually the 
efforts focused on the understanding of the disturbing ef- 
fects of the finite length of the sequences on the estimate 
of the fractal dimension Dt as well as on the multifractal 
characterization of the trail. For instance, although it 
is well-known that in two-dimensions an infinite length 
true random walk (i.e., a sequence of bases generated at 
random) is space filling and so Dt = 2, the box-counting 
and sandbox methods yield estimates significantly lower 
than 2 even for walks as large as 10^ steps ||l^. Moreover, 
finite-length random walks exhibit an effective multifrac- 
tal spectrum due mainly to crossover effects between the 
usually distinct bulk and surface properties of the trail 
Ip^ . These caveats, however, have not prevented claims 
that relatively short (typically 2 x 10'' bp) mitochondrial 
DNA genomes have a definite multifractal structure . 
In Sec. Ill of this contribution we invoke a classical result 



of fractal theory to link the Hurst exponent of the DNA 
walk to the fractal dimension of the corresponding trail 
in space (see, e.g., |o|), then arguing in favor of the 
monofractality of DNA trails too. 

II. DNA WALKS 

Consider a record, such as that given in Eq. (|l|), where 
the variables yi are not necessarily Ising variables but 
must form a compact support, i.e., yi = for a very 
small fraction of the elements only. The generalized mul- 
tifractal DFA (MF-DFA) involves the following steps (we 
refer the reader to ref. for a thorough presentation 
of the method). First, divide the entire sequence into 
Ni = v[it{N/l) nonoverlapping segments of length /. Sec- 
ond, for each segment v = 1, . . . ,Ni of length I calculate 
the local trend by a least-square fit of the record in the 
segment, and denote by Si,{n) the fitting polynomial in 
segment v. Then evaluate the variance 

^) = y E [('^ - 1) ^ + - ^''H}' (4) 

71=1 

for each segment. An important parameter is the order 
TO = 0, 1, 2, ... of the polynomial Sj^(n) used in the fitting 
procedure. The choices of to correspond to different or- 
ders of DFA, denoted by DFAm, which differ in their ca- 
pability of eliminating trends in the series. For instance, 
DFAO can eliminate only constant trends in the series, 
DFAl eliminates constant as well as linear trends, and 
so on, so that a comparison of the results for different 
orders of DFA yields information on the type of trend 
in the series. Third, determine the order q fluctuation 
function averaging over all segments 

where in general the index variable q takes on any real 
value, except zero. As the final step, analyze log-log plots 
of Fq{l) versus I for each value of q in order to determine 
the scaling behavior of the fluctuation functions. For 
large values of I, we expect these functions to increase 
with increasing I as a power law 

F,il)^l''('^l (6) 

Clearly, for q = 2 one has h{2) — H hy construction, and 
in the case that the scaling behavior of F'^{1, v) is iden- 
tical for all segments i.e., the record is a monofractal, 
the exponent hiiq) is independent of as expected. The 
definition (||) is well-suited to detect discrepancies in the 
scaling behavior of the large and small fluctuations. In 
particular, for positive g, ^.(g) yields the scaling behavior 
of the segments with large fluctuations while for negative 
g, /i(g) yields the scaling of the segments with small vari- 
ances. In addition, there is a simple relation between h((i) 




q 

FIG. 3: Multifractal spectrum r(g) vs q for the DNA 
walk of the intron-rich HUMHBB sequence depicted in Fig. 
|l[ The data correspond to different orders of DFA, m = 
(0)i 1 (V)i 3 (x) and 5 (A). The symbol + corresponds to 
the data for the shuffled sequence. 

and the scaling exponents T(g) defined by the standard 
partition function-based multifractal formalism [ p^ , 

T{q) - qh{q) - 1. (7) 

The importance of this relation is that it allows compar- 
ison of the results obtained in the MF-DFA scheme with 
those of the standard multifractal analysis in the case of 
stationary series, and with those of the WTMM in the 
case of nonstationary series. 

In Fig. 1^ we show the multifractal spectrum r(g) ob- 
tained by the application of different orders of the MF- 
DFA method to the intron-rich HUMHBB sequence. The 
results for orders larger than 5 are indistinguishable in 
the scale of the figure, so the data for m = 5 (A) yields 
the trendless spectrum. Analysis of this figure gives valu- 
able information on the nature of the correlations of the 
DNA walk. In particular, the linearity of the trendless 
spectrum indicates that the record s(n) depicted in Fig. 
|l| is monofractal. Actually, this conclusion is not affected 
by the presence of trends in the series, since the q depen- 
dence of T(g) is linear regardless of the value of m. How- 
ever, the value of the scaling exponent H = h{q) Vg, given 
by the slope of the straight lines that fit the data, is sensi- 
tive to local trends. For example, we find H — 0.72±0.02 
for TO = 0, i7 = 0.68±0.02 for to = 1, and = 0.60±0.02 
for TO = 5. As expected, H = 0.50 ± 0.02 for the shuf- 
fled sequence. All these numerical values are in perfect 
agreement with the values obtained using the WTMM 
method Q . Since the trendless scaling exponent is defini- 
tively different from 1/2 we can conclude for the exis- 
tence of long-range correlations in this particular intron- 
containing sequence. The long-range correlation induced 
by the local trends is reflected in the larger value of the 
scaling exponent calculated with the MF-DFAO, which 
overestimates the value of H in about 20%. It is in- 
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FIG. 4: Same as Fig. | but for the DNA walk of the E. 
coli chromosomal sequence ECOllOK composed primarily of 
coding regions. 



out using the WTMM method, for several representa- 
tions proposed in the hterature (see, e.g. |^, |l^ ^M)- 
As we will show in the sequel, the results were essentially 
the same as those reported here for the one-dimensional 
DNA walks. 



III. DNA TRAILS 

As pointed out before, a DNA trail is the projection of 
a two-dimensional DNA walk into the space plane. Fig. 
H shows the time record of the walk that gives rise to the 
trail illustrated in Fig. ||for the HUMHBB sequence. We 
begin by noting that the generalization of the MF-DFA 
method to two or higher-dimensional records is straight- 
forward as one needs only to change the basic Eq. (jj), 
which is rewritten as 



1 ' 

^) = y ^ {s - 1) Z + n] - s„(n)}2 (8) 



teresting to note that the trends influencing the large 
fluctuations {q > 0) are practically unaffected by the ap- 
plication of DFAl. 

The results of the application of the MF-DFA method 
to the characterization of a sequence composed predomi- 
nantly of coding regions is summarized in Fig. ^ The se- 
quence considered is a portion of the E. Coli K12 genome 
(GenBank name ECOllOK) for which N = 111401. We 
find H = 0.60 ± 0.02 for m = 0, and iJ = 0.51 ± 0.02 for 
m > 1 as well as for the shuffled sequence. As before, the 
monofractal structure of the record is confirmed by the 
linearity of the T{q) spectrum. Moreover, failure to elim- 
inate the local trends leads to a wrong conclusion about 
the existence of long-range correlations in this DNA walk, 
as evidenced by the overestimate of H that results from 
the application of MF-DFAO. The local trends in this se- 
quence seem to have a particularly simple linear nature 
since application of MF-DFAl is sufficient to eliminate 
them altogether. 

To conclude this section we note that we have consid- 
ered several alternative one-dimensional representations, 
such as the strong-weak bond classification where C or 
G are associated to one step up {yi = 1) and A or T to 
one step down {yi — —1). All the alternative schemes 
analysed have produced within the numerical precision 
exactly the same exponents of the purine-pyrimidine clas- 
sification. Our results could also be obtained within the 
WTMM framework and, in fact, the wavelet analysis of 
the HUMHBB sequence can be found in ref. . The aim 
of this section was to give additional support to the claim 
of Kantelhardt et al. that, despite the operational and 
conceptual simplicity of the MF-DFA method, it yields 
results that are identical to those derived by the sophis- 
ticated WTMM method |l^. For instance, this easi- 
ness of implementation has allowed us to study higher- 
dimensional DNA walks, an awkward task to be carried 



where s(n) = x(n)i + y(n)j and similarly for the fit- 
ting polynomial vector s^. Over an interval of I steps 
the position vector in the trail will vary by typically 
I 6r |« F2{1) = I" while the 'mass' M of the trail, i.e., 
the number of points generated by these steps is M ^ I, 
provided the contribution of overlapping points is negli- 
gible. Hence M ^\ Sr |^/^, from where one concludes 
that Dt — 1/H. This relation is valid only if l/_ff < d 
where d is the space dimension. When l/-ff > d overlap 
cannot be neglected and so Dt = d (see 1^, |l^ for math- 
ematical and pictorial details of this argument). Thus 
knowledge of the Hurst exponent H for the d-dimcnsional 
DNA walks determines uniquely the fractal dimension of 
the DNA trail. Calculation of H or, more generally, of 
h{q) follows the same procedure sketched before for the 
one-dimensional case, except that the variance given in 
Eq. (^ is now replaced by F^. 

Appfication of the MF-DFA method to two- 
dimensional DNA walks representing the HUMHBB (see 
Fig. |) and the ECOllOK sequences yields, within the 
numerical precision, the same results as for the one- 
dimensional case. Actually, so far as the value of the 
Hurst exponent H is concerned this agreement is ex- 
pected since using Eqs. (|) and (|) for g = z we 
have simply F2 ^ l^^ + l^" , where and Hy are the 
scaling exponents of the one-dimensional DNA walks in 
the plane {n,x) and (n,y). Hence for large I one has 
H — max{Hx, Hy). In Table | wc present the estimate of 
the fractal dimension Dt of the HUMHBB and ECOllOK 
trails, as well as of a random trail characterized by the 
same length and base frequency as the ECOllOK se- 
quence. The errors in the estimates of Dt are statistical 
errors. For the box-counting method the systematic er- 
rors are probably much larger than those shown in Table 
|, due mainly to finite- length effects [ p^ . 

Interestingly, although the Dt estimates of the MF- 
DFArn method are based on the analysis of finite-length 
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FIG. 5: Two-dimensional DNA-walk representation of the 
HUMHBB sequence whose projection into the space plane 
yields the trail of Fig. |2|. 



Sequence 


N 


MF-DFA5 


MF-DFAO 


box-counting 


HUMHBB 
ECOllOK 
random 


73326 
111401 
111401 


1.67 ±0.06 
1.96 ±0.08 
2.00 ±0.08 


1.39 ±0.04 
1.67 ±0.06 
2.00 ±0.08 


1.40 ±0.01 
1.38 ±0.01 
1.45 ±0.01 



TABLE I: Estimates of the fractal dimension Dt of two- 
dimensional trails. The MF-DFAm estimates are based on 
the formula Dt = ^IH where H is the Hurst exponent of the 
record (e.g. Fig. H). The box-counting method is applied 
directly to the trail (e.g. Fig. 



records, this has practicaUy no effect on the results, as ev- 
idenced by the correct calculation of Dt for the random 
sequence, a test the box-counting method clearly fails. 
Moreover, these data corroborate the intuition, stemming 
from the visual inspection of Fig. |[ that the effect of the 
local trends is to decrease the fractal dimension of the 
trail. The excellent agreement between the box-counting 
estimate, which is affected by both local trends and finite- 
length effects, and the MF-DFAO estimate, which is af- 
fected by local trends only, indicates that the surface ef- 
fect is negligible for the HUMHBB sequence. In fact, this 
is expected since the box-counting method is very reliable 
to calculate the fractal dimensions of structures that do 
indeed exhibit fractality, predictably failing in the case 
of space-filling structures. The results for the ECOllOK 
sequence illustrate the bad performance of box-counting 
when the two sources of systematic errors are present. 
The aforementioned agreement between the box-counting 
and the MF-DFAO methods demonstrates inequivocally 
that the spurious effects of local trends must be filtered 
out for a correct account of the fractal or multifractal 
properties of DNA trails, a procedure that has been ig- 
nored in all previous analyses presented in the literature 



Since the MF-DFA analysis has indicated that even 
the higher-dimensional DNA walks are monofractals, one 
concludes that the variance F"^ {l,^) defined in Eq. (||) 
is the same for all segments v of the record. As a re- 
sult, the variation of the position vector | | in the 
trail will indcpend of the segment v and so the trail must 
have a monofractal structure too. This conclusion seems 
at odds with those of multifractal analyses based on the 
sandbox and box-counting methods ^|. How- 

ever, we note that the original claim of ref. |l6| was 
that finite-length DNA trails, as well as randomly gen- 
erated trails, are characterized by effective multifractal 
spectra. Clearly, these spectra are artifacts of the sand- 
box and box-counting methods which cannot deal ade- 
quately with the surface effects. The unjustified dropping 
of the adjective 'effective' ^S), led then to the appar- 
ent disagreement between the two approaches. These 
effective multifractal spectra might indeed be an useful 
tool for comparing random or biological sequences of sim- 
ilar lengths ijl^, especially if one manages to free them 
from the local trends effects. 



IV. CONCLUSION 

The main purpose of this contribution was to point 
out that the statistical properties of two apparently dis- 
tinct and pictorially appealing representations of DNA 
sequences, namely, DNA walks (Figs. |l| and ||) and DNA 
trails (Fig. ||), are in fact closely related. In particular, 
we argued that the patchiness of the DNA sequences, 
which has greatly hindered the characterization of the 
long-range correlations of DNA walks || , also plagues 
the fractal and multifractal analyses of the DNA trails, a 
fact that has been largely overlooked in previous investi- 
gations |l^, Using the recently proposed multi- 
fractal detrended fluctuation analysis (MF-DFA) method 
]l3[ to filter the local trends of DNA walks in d = 1 and 
2 dimensions, we were able to show that these walks are 
monofractals. This conclusion holds true for the non- 
detrended walks as well, since the multifractal spectrum 
T{q) is linear regardless of the order m > of the fitting 
polynomials used to eliminate the local trends (Figs. ^ 
and|^). More importantly, the independence of the vari- 
ance F'^{l,v) [see Eqs. (^) and (^)] on the segment v 
of the DNA walk carries over to the DNA trail which 
is thus monofractal too. In particular, the fractal di- 
mension of the trail is Dt = 1/ H where H is the Hurst 
exponent of the record 10 . In the case 1/H > d, 
as for the one-dimensional walks, one has Dt ~ d. Al- 
though we have presented data for the HUMHBB and 
the ECOllOK sequences only, which probably are, re- 
spectively, the most popular examples of intron-rich and 
intron-poor sequences, we have verified that our main re- 
sults hold true for many other sequences in the GenBank. 

Once one has discarded the simple local variations 
in nucleotide composition along DNA sequences as the 
cause of the long-range correlations observed in intron- 



6 



rich sequences (though these variations do cause the spu- 
rious long-range correlations in intron-poor sequences), 
it is natural to ask then what are the sources for these 
correlations, which, as shown above, are also respon- 
sible for the nontrivial (fractal) geometry of the DNA 
trails. The answer is that neither the internal structure 
of patches nor their order in the sequence are relevant: it 
is the power-law distribution of patch lengths that deter- 
mines the true long-range correlations Q. This as well 
as several other findings concerning the statistical prop- 
erties of DNA sequences have prompted the proposal of 
minimal evolutionary models, based on biologically mo- 
tivated mechanisms, to account for those features (see, 
e.g., ||2^, |2^, HI). In fact, rather than providing mea- 
sures to characterize or distinguish classes of sequences, 
we think the thrust of the research on large-scale sta- 



tistical properties of genomes is to provide quantitative 
standards for modeling the inherently stochastic process 
of molecular evolution p^ , bearing thus on fundamen- 
tal issues such as the origin of life and the evolution of 
complexity p4| . 



Acknowledgments 

This research was supported by Fundagao de Amparo 
a Pesquisa do Estado de Sao Paulo (FAPESP), Proj. No. 
99/09644-9. The work of J.F.F. was supported in part 
by Conselho Nacional de Desenvolvimcnto Cientffico e 
Tecnologico (CNPq). 



[1] C.-K. Peng, S. V. Buldyrev, A. L. Goldberger, S. Havlin, 
F. Sciortino, M. Simons, and H. E. Stanley, Nature 356, 
168 (1992). 

[2] C.-K. Peng, S. V. Buldyrev, S. Havlin, M. Simons, H. E. 
Stanley, and A. L. Goldberger, Phys. Rev. E 49, 1685 
(1994). 

[3] S. Karlin and V. Brendel, Science 259, 677 (1993). 

[4] C. K. Chui, An Introduction to Wavelets (Academic 
Press, Boston, 1992). 

[5] P. Bernaola-Galvan, R. Roman-Roldan, and J. L. Oliver, 
Phys. Rev. E 53, 5181 (1996). 

[6] A. Arneodo, E. Bacry, P. V. Graves, and J. F. Muzy, 
Phys. Rev. Lett. 74, 3293 (1995). 

[7] A. A. Tsonis, P. Kumar, J. B. Eisner, and P. A. Tsonis, 
Phys. Rev. E 53, 1828 (1996). 

[8] L. Luo, W. Lee, L. Jia, F. Ji, and L. Tsai, Phys. Rev. E 
58, 861 (1998). 

[9] B. B. Mandelbrot, The Fractal Geometry of Nature (Free- 
man, San Francisco, 1982). 
[10] R. F. Voss, Physica D 38, 362 (1989). 
[11] J. Feder, Fractals (Plenum Press, New York, 1988). 
[12] T. C. Halsey, M. H. Jensen, L. P. Kadanoff, I. Procaccia, 



and B. L Shraiman, Phys. Rev. A 33, 1141 (1986). 
[13] J. W. Kantelhardt, S. A. Zschiegner, E. Koscielny- 

Bunde, A. Bunde, S . Havlin, and H. E. Stanley, 

arXiv:physics/0202070 , 
[14] M. A. Gates, J. theor. Biol. 119, 319 (1986). 
[15] C. L. Berthelsen, J. A. Glazier and M. H. Skolnick, Phys. 

Rev. A 45, 8902 (1992). 
[16] C. L. Berthelsen, J. A. Glazier, and S. Raghavachari, 

Phys. Rev. E 49, 1860 (1994). 
[17] J. A. Glazier, S. Raghavachari, C. L. Berthelsen and M. 

H. Skolnick, Phys. Rev. E 51, 2665 (1995). 
[18] N. N. Oiwa and J. A. Glazier, Physica A 311, 221 (2002). 
[19] T. Vicsek, Fractal Growth Phenomena (World Scientific, 

Singapore, 1992). 
[20] W. Li, Phys. Rev. A 43, 5240 (1991). 
[21] P. M. C. de Oliveira, Physica A 273, 70 (1999). 
[22] Y. Almirantis and A. Provata, BioEssays 23, 647 (2001). 
[23] M. Kimura, Theor. Pop. Biol. 2, 174 (1971). 
[24] J. T. Bonner, The Evolution of Complexity (Princeton 

University Press, Princeton, 1988). 



